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This  paper  summarizes  the  results  presented  at  the  Army  Research  Workshop  held 
at  Monterey,  CA  in  October,  1989.  A  more  detailed  version  will  appear  elsewhere. 

—  ^  In  the  age-dependent  minimal  repair  model  of  Block,  Borges,  and  Savits  (1985), 
a  system  failing  at  age  t  undergoes  one  of  two  types  of  repair.  With  probability 
p(t),  a  perfect  repair  is  performed,  and  the  system  is  returned  to  the  *good-as-neWr' 
state,  while  with  probability  1  -  p(t),  a  minimal  repair  is  performed,  and  the  sys¬ 
tem  is  repaired,  but  is  only  as  good  as  a  working  system  of  age  t.  Whitaker  anc 
Samaniego  (1989)  propose  an  estimator  for  the  system  life  distribution  F  when  data\iNSPtcttD^ 
are  collected  under  this  model. 

Using  the  product  integral  representation  of  the  survival  function,  a  basic  result 
of  Block,  Borges,  and  Savits  concerning  the  waiting  time  until  the  first  perfect  repair 
is  extended  to  allow  for  discontinuous  distribe  Jons.  Then  using  counting  process 
techniques,  the  large  sample  theorems  of  Whitaker  and  Samaniego  are  extended  to 
the  whole  line.  These  results  are  used  to  derive  confidence  bands  for  F,  and  to 
determine  a  sufficient  condition  for  their  applicability  on  the  whole  line.  Simulation 
results  for  the  bands  are  provided.  An  extension  of  the  Wilcoxon  two-sample  test  to 
the  minimal  repair  model  is  also  examined.  !  i  r- 


1  The  Minimal  Repair  Model 


To  fix  notation,  let  F  be  a  life  distribution,  let  tj?  be  the  upper  endpoint  of  the  support  of  F 
(possibly  infinite),  and  let  A(f)  =  /(o,<](F(s  —  ))~1dF(s)  be  the  cumulative  hazard  function 
of  F,  where  F  =  1  —  F. 

Now,  for  j  =  1, . . .  ,n,  let  {  Xjto  =  0,  Aj,i ,  Xj<2,  •  •  •}  be  independent  record  value 
processes  from  F.  These  are  Markov  processes  with  P{Xj t*  >  t  j  Xj$, . . . ,  A' y,fc_i )  = 
F(t)/F(Xj'k~i) ,  for  t  >  Xjtk- 1,  k  >  1.  If  A F(tf)  >  0,  define  Xjj  =  co  for  all  /  larger  than 
the  first  k  for  which  Xjtk  =  Tp.  In  all  cases  we  take  p(rp )  =  1.  These  processes  represent 
the  failure  ages  of  n  systems  under  a  “forever  minimal  repair”  scheme. 
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Perfect  repair  is  introduced  into  this  model  by  the  use  of  independent  uniform  random 
variables.  This  facilitates  the  construction  of  the  a-field  structure  (filtrations)  necessary 
to  our  analysis  of  the  model  through  martingale  methods.  Thus  we  let  {  Uj,k  ■  1  <  j  < 
n,  k  >  1  },  be  i.i.d.  uniform  r.v.’s,  and  define 

Sj'k  =  I(Ujtk  <  p(Xjtk)), 

Vj  =  inf{  k  :  Sjik  =  1  }. 

Thus  observing  { (Ayi, . . . ,Xj^3);j  =  l,...,n),  is  equivalent  to  observing  n  indepen¬ 
dent  copies  of  the  age-dependent  minmal  repair  process  of  Block,  Borges,  and  Savits 
(BBS)(1985),  each  until  the  time  of  its  first  perfect  repair. 

This  structure  provides  us  with  a  concrete  starting  point  for  a  statistical  analysis  of  the 
BBS  model.  However,  we  need  conditions  which  are  sufficient  to  assure  the  finiteness  of 
Xj,u}  •  Such  conditions  are  given  by  the  following  result,  which  generalizes  a  result  of  BBS  to 
the  case  of  possibly  discontinuous  F.  Though  this  generalization  may  not  be  important  for 
modeling  system  failures,  it  will  be  useful  to  us  in  proving  large  sample  results.  Also,  the 
proof  of  this  result,  which  we  sketch  below,  is  more  straightforward  than  the  original  proof 
of  BBS.  The  reader  ’S  referred  to  Hollander,  Proschan,  and  Sethuraman  (1989)  (HPS),  for 
detailed  proofs  of  this  and  other  results  in  this  paper. 

Proposition  1  Let  H(t)  =  P(XV  <t,v<  oo).  Then 

m  =  n(o,«](i  -  dA„) 


Moreover,  if  either 

(i)  AF(tf)  >  0  ( and  p(rF)  =  1), 

or 

(ii)  F(tf- )  =  1  and  =  +00, 

then  H  is  a  proper  distribution  function  and  v  is  almost  surely  finite.  Conversely,  if  H  is 
a  proper  distribution  function,  then  either  (i)  or  (ii)  must  hold. 

Proof.  (Sketch)  Note  that 

H(t)  =  1  -P(X„  <  t,u<  oo) 

=  i -Em,  <<,*'=.>)• 

i=i 

A  conditioning  argument  shows  that 

+  f  /  ••/  Mh)--  da(l y), 
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where 


This  is  equivalent  to 


a (0  =  /  (1  -Pi*)) 


dF(s) 
F(s)  ■ 


pttt  =  iw1  + da )  =  exp  (aC(0)  nc1  +  Aq(5)). 

r  VV  3<t 

where  ac  is  the  continuous  part  of  a  and  Aa(t)  is  the  jump  in  a  at  t.  Here,  E[(o,(](l  +  da) 
represents  a  product  integral.  The  theory  of  product  integration  with  applications  in 
statistics  is  reviewed  in  Gill  and  Johansen  (1987).  The  result  follows  from  the  last  equation 
after  some  algebra.  □ 

We  will  say  that  a  pair  satisfying  either  (i)  or  (ii)  describes  a  regular  repair  scheme. 


2  The  Whitaker- Samaniego  Estimator 

In  this  section,  we  derive  a  martingale  representation  for  the  Whitaker-Samaniego  (1989) 
estimator  (WSE).  This  representation  is  then  used  in  conjunction  with  Rebolledo’s  Mar¬ 
tingale  Central  Limit  Theorem  and  the  techniques  of  Gill  (1983)  to  derive  limit  theorems 
for  the  WSE. 


The  Basic  Martingale 

Define 

N*(t)  =  #{k:Xj<k<t}, 

and 

Ft  =  v  ({  NJ(s)  :  s  <  t,  1  <  j  <  n  }) 

V  <r({Ujtk:k>  1,1  <j  <n}). 


For  the  rest  of  this  paper,  (^)t>o  will  serve  as  the  underlying  filtration  for  all  martingales. 
Now  let 

N(t)  =  #{  ( j ,  k)  :  Xj'k  <t,k<  1  <  j  <  n  }, 

Y(t)  =  #{;'  ••  X} -tVj  >  t,  1  5;  i  :  n  }, 

and 

M(t)  =  N(t)~ j  Y(s)d\(s). 

In  HPS,  it  is  shown  that  M  is  a  locally  square-integrable  martingale  with  predictable 
quadratic  variation  given  by 

(M)(i)=  f  Y(s)(l-AA(s))d\(s).  (1) 


This  provides  the  basic  martingale  structure  for  further  analysis  of  the  minimal  repair 
model. 


A  Martingale  Representation  for  the  WSB 

Assume  that  F  is  continuous  and  that  the  pair  ( F,p )  describes  a  regular  repair  scheme. 
Let  X(k)  be  the  kfi1  ordered  value  of  the  set  {  Xjtk  '■  k  <  vj,  1  <  j  <  n  },  let 

T  =  min{  X(k)  :  Y(X(k))  =  1  }, 


and  let  J(s)  =  I(s  <  T).  Then  the  Whitaker-Samaniego  estimator  (WSE)  can  be  written 
as 

m = n(o.,)(i  -  iK)  =  n  (i  -  aA(»>)  , 

3<t 

where 

Using  Duhammel’s  equation  (Gill  and  Johansen,  1989),  (F  —  F)/F  can  be  expressed 
as  an  integral  with  respect  to  the  martingale  M: 


m  -  m 
m 


t*  h*-) 

Jo  F(s)Y(s) 


dM(s). 


From  this  and  (1)  it  follows  (F  —  F)/F  is  itself  a  locally  square-integrable  martingale  with 
predictable  quadratic  variation  process  given  by 

/LzZ\  -  [  (h>-)y  m*) 

\  P  /  Jm\F(s-)J  F(s)Y(s)' 

This  quadratic  variation  process  essentially  serves  to  identify  the  covariance  structure  of 
the  limiting  Gaussian  processes  derived  in  the  next  section. 


Large  Sample  Results 

With  the  above  representation,  Rebolledo’s  martingale  CLT  and  the  methods  of  Gill(  1 983) 
yield  the  following  result,  which  extends  Theorem  3.3  of  Whitaker  and  Samaniego  (1989) 
to  the  whole  line. 

Theorem  1  Let  ( F,p )  describe  a  regular  repair  scheme,  with  F  continuous.  Then  the 
following  hold: 

(i)  As  n  — >  oo, 

y/n(F  -  F)  £  F-B{C)  in  Z>[0,oo], 
where  B  is  Brownian  motion  on  [0,oo),  and 

ca )  =  r  dF m 

U  Jo  H(s-)F{s)' 


4 


(ii)  As  n  -»  oo, 

y/^^(F  -  F)  Z  B°{I<)  m£>[0,oo], 

F 

where  B°  is  Brownian  bridge  on  [0,1],  and  K  =  C/(  1  +  C). 

Details  of  the  proof  of  this  theorem  are  given  in  HPS.  We  note  here  that  the  proof  of 
(i)  does  not  require  any  additional  conditions  beyond  regularity  of  the  repair  scheme.  This 
is  in  contrast  with  the  analogous  result  of  Gill  (1983)  for  the  Kaplan-Meier  estimator  in 
the  usual  censored  survival  data  model,  where  some  condition  on  the  amount  of  censoring 
is  needed.  We  will  see  below  however,  that  an  additional  condition  limiting  the  amount  of 
imperfect  repair  is  needed  to  assure  convergence  of  the  expression  in  (ii)  when  an  estimate 
is  substituted  for  K/F. 


3  Applications 

In  this  section,  the  asymptotic  results  of  the  last  section  are  used  to  derive  large  sample 
confidence  bands  for  F  and  to  obtain  the  limiting  distribution  of  an  extension  of  the 
Mann-Whitney-Wilcoxon  test  statistic  to  the  minimal  repair  model. 


Confidence  Bands 

The  result  in  part  (ii)  of  Theorem  1  suggests  confidence  bands  based  on  the  distribution 
of  the  supremum  of  Brownian  bridge.  It  is  necessary  however  to  estimate  K/F  in  order 
to  construct  the  bands.  Let  H  be  the  empirical  cdf  of  the  XjyU] ,  and  let  I<  =  C/(  1  -f  C), 
where  C  is  defined  by 

<5(0  =  /  . iF(>l  , 

Am  h{s-)F(s ) 

We  would  like  to  have 

y/n  r~-(F  —  F)  — ►  B°(K)  in  D[0,  oo],  as  n  — ►  oo,  (2) 

F 

in  order  to  justify  asymptotic  (1  —  a)  x  100%  confidence  bands  for  F  of  the  form 

F±y/H\aF/K, 

where  \a  is  the  upper  quantile  of  the  distribution  of  sup  |£°(f)|. 

We  can  show  that  (2)  holds  on  [0,  r]  for  any  r  <  Tp,  but  for  the  complete  result,  some 
additional  condition  seems  to  be  needed.  Using  the  result  of  Prop.l,  it  is  shown  in  HPS 
that  K/F  and  K/F  are  nondecreasing  and  that 
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and 
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Using  this,  it  can  be  shown  that  a  sufficient  condition  for  (2)  is  that 


\=[  (1-P(s))d\{s)  <oo. 

F(tf-)  J(  o,tf) 

This  condition  requires  that  p(t)  — >  1  as  t  |  tf  (at  a  rate  sufficient  for  the  convergence  of 
the  integral),  and  hence  provides  a  limit  on  the  amount  of  imperfect  repair. 

Simulation  results  for  the  bands  computed  over  finite  intervals  (in  the  case  of  constant 
p )  indicate  that  coverage  probabilities  are  quite  good  for  sample  sizes  of  50  or  more.  This 
will  of  course  vary  with  the  parameters  of  the  model.  Simulations  were  carried  out  with 
both  Gamma  and  Weibull  F ,  with  varying  shape  parameters,  and  with  various  values 
of  p,  various  interval  lengths,  and  various  nominal  confidence  levels.  As  an  example,  the 
following  table  gives  the  simulated  coverage  probabilities  for  nominal  95%  confidence  bands 
over  the  interval  [0, 4.744]  when  the  underlying  F  is  Gamma  with  shape  parameter  2.  (Note 
that  4.774  is  the  ninety-fifth  percentile  of  Gamma(2).)  More  extensive  tables  are  provided 
in  HPS. 


n 

p  =  .50 

p  =  .25 

p  =  .10 

10 

.9025 

.8660 

.8710 

20 

.9270 

.9125 

.9187 

30 

.9460 

.9287 

.9327 

50 

.9515 

.9398 

.9395 

100 

.9528 

.9540 

.9452 

200 

.9515 

.9517 

.9495 

An  Extension  of  the  Mann-Whitney-Wilcoxon  Test 

Using  part  (i)  of  Theorem  1,  it  is  also  possible  to  obtain  the  limiting  distribution  for 
an  adaptation  of  the  Mann-Whitney-Wilcoxon  two-sample  statistic  to  the  minimal  repair 
model.  Here  we  assume  that  for  i  =  1,2,  we  observe  nj  BBS  processes  from  (F,-,p,),  each 
until  its  first  perfect  repair.  In  general  we  wish  to  test  the  null  hypothesis  H0 :  F\  =  F2, 
with  typical  one-sided  alternatives  specifying  JFidF2  >  1/2,  and  two-sided  alternatives 
specifying  / dF \  ^  1/2. 

A  statistic  analogous  to  the  Mann- Whitney  form  of  the  Wilcoxon  two-sample  test 
statistic  is  VU,  as  given  by 


W 


E 

<w2(*)>o 


Fi(i >)F2{s~) 


A  N2(s) 
Y7(S)  ’ 


where  F{  is  the  WSE,  A JV,(s)  is  the  number  of  failures  at  age  s,  and  T(s)  is  the  num¬ 
ber  of  items  at  risk  at  age  s  in  the  ith  sample.  This  statistic  is  a  natural  estimator 
of  JF-[dF 2  =  P(X i  <  X2),  where  X\  and  X2  are  independent  random  variables,  with 
Xi  ~  F{.  Assuming  continuous  distributions,  P{X\  <  A^)  =  1/2  under  Ho,  and  in  the 
one-sided  case,  significantly  large  values  of  W  provide  evidence  against  Ho  in  the  direction 
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of  f  Fi  (IF2  >  1/2.  For  large  sample  sizes,  we  have  the  following  result,  which  is  proven  in 
HPS: 


Theorem  2  If  Fi  and  F2  are  continuous,  and  the  pairs  (F\,pi)  and  (F2,P2)  describe  reg¬ 
ular  repair  schemes,  and  if  n\,ri2  — ►  00  in  such  a  way  that  ninfn2  — ►  A,  0  <  A  <  1, 
then 

r  r  1 _  /  1  1  \ 

(3) 


where 


vATTF ~2^W-fF1dF2]  2  N(o,  + -Ljaiy 

noo 

F1(s)F1(t)C1(t)dF2(s)dF2(t), 

noo 

F2(s)F2(t)C2(t)dF1(s)dF1(t). 


Under  the  null  hypothesis,  Ho  :  F\  =  F  =  F2 , 


=  2  £  F(t)c,(t)  (£  F(s)JF(s))  dF(t )  =  i  £  rj£L  <if{s > . 

For  purposes  of  testing  the  null  hypothesis  in  the  large  sample  case,  we  thus  propose 
referring  the  test  statistic 


2' 

2 

n2 


to  a  standard  normal  distribution,  where 


1  r  £d?Ldp.(S)  =  xl  £ 

Ht(s~)  4 


nF  (s)Fi(s-)ANi(s) 
Y?(s) 


ANi(s)X) 

and  Hi  is  the  empirical  distribution  of  the  perfect  repair  ages  in  the  ith  sample. 

It  is  shown  in  HPS  that  the  a,  are  consistent,  which  justifies  the  use  of  this  test.  If 
the  pi  are  constants  (see  Brown-Proschan  (1983)),  the  above  expressions  simplify  greatly 
under  H0.  If  Fi  =  F2  =  F,  then  Hi  =  FPi ,  and  the  asymptotic  variance  in  (3)  reduces  to 

1  2  1  2  if  1  \  t  1  (  1  ^ 

+  1  _  A"2  "  A  U*  -  Pi) )  1  -  A  \4(4  -  p,) )  ' 

The  pfs  are  of  course  consistently  estimated  by  their  MLE’s,  p,-,  the  ratio  of  n,  to  the  total 
number  of  failures  in  the  zth  sample,  and  for  large  samples,  the  statistic  Z',  given  by 

1/2 


z'=(w-  5)/ 


1 


+ 


1 


[4n1(4  —  Pi)  4n2(4-p2). 

can  be  referred  to  a  standard  normal  distribution  in  order  to  test  the  null  hypothesis.  Note 
also  that  if  pi  =  p2  =  1,  then  we  are  in  the  usual  i.i.d.  two-sample  model,  the  WSE’s 
reduce  to  the  empirical  c.d.f.’s,  and  W  is  just  a  multiple  of  the  Mann- Whitney  form  of  the 
Wilcoxon  rank-sum  statistic.  In  this  case,  the  above  results  yield 

in  agreement  with  the  usual  results  for  the  Mann-Whitney-Wilcoxon  test. 
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This  paper  summarizes  the  results  presented  at  the  Army  Research 
Workshop  held  at  Monterey,  CA  in  October,  1989.  A  more  detailed  version 
will  appear  elsewhere. 

In  the  age-dependent  minimal  repair  model  of  Block,  Borges,  and 
Savits(1985) ,  A  system  failing  at  age  t  undergoes  one  of  two  types  of  repair. 
With  probability  p(t),  a  perfect  repair  is  performed,  and  the  system  is 
returned  to  the  "good-as-new"  state,  while  with  probability  l-p(t),  a  minimal 
repair  is  performed,  and  the  system  is  repaired,  but  is  only  as  good  as  (ctd) 
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20.  (ctd)  a  working  system  of  age  t.  Whitaker  and  Samaniego  (1989)  propose  an 
estimator  for  the  system  life  distribution  F  when  data  are  collected  under  this 
model . 

Using  the  product  integral  representation  of  the  s  rvival  function,  a  basic 
result  cf  Block,  Borges,  and  Savits  concerning  the  waiting  time  until  the  first 
perfect  repair  is  extended  to  allow  for  discontinuous  distributions.  Then, 
using  counting  process  techniques,  the  large  sample  theorems  of  Whitaker  and 
Samaniego  are  extended  to  the  whole  line.  These  results  are  used  to  derive 
confidence  bands  for  F,  and  to  determine  a  sufficient  condition  for  their 
applicability  on  the  whole  line.  Simulation  results  for  the  bands  are  provided. 
An  extnesion  of  the  Wilcoxon  two-sample  test  to  the  minimal  repair  model  is  also 
examined. 


